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di! Abstract 

' New, radiatively generated, LO quark (u, d, s, c, b) and gluon densities in a real, unpolarized 

^ - photon are presented. We perform a global 3-parameter fit, based on LO DGLAP evolution 

. equations, to all available data for the structure function F2{x,Q'^). We adopt a new theoretical 

approach called ACOT(x), originally introduced for the proton, to deal with the heavy-quark 
thresholds. This defines our basic model (CJKL model), which gives a very good description of 
the experimental data on for both and x dependences. For comparison we perform 

a standard fit using the Fixed Flavour-Number Scheme (FFNSc ji^L model), updated with respect 
to the previous fits of this type. We show the superiority of the CJKL fit over the FFNSc one 
and other LO fits to the i^2^(x,(5^) data. The CJKL model gives also the best description of the 
LEP data on the dependence of the F2 , averaged over various x-regions, and the Fg"^, which 
were not used directly in the fit. Finally, a simple analytic parametrization of the resulting parton 
densities obtained with the CJKL model is given. 
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I. INTRODUCTION 



The photon structure function was recognized as an important quantity aheady in the 
early days of the Parton Model It has attracted even more attention since the seminal 
paper by Witten ^, which shows that F2 can serve as a unique test of QCD. This expec- 
tation was based on the fact that the (asymptotic) point-like solution of the evolution 
equation, summing the leading QCD corrections, can be obtained for F2 without additional 
assumptions. Further studies showed the need of the hadronic, VMD-type, contribution to 
F^, and consequently the need of an input, as for every other hadronic structure function. 

The structure function F2 is extracted from measurements in deep inelastic scattering 
on a photon target, which can be performed in e^e~ experiments. The data are used 
to construct parametrizations of the parton distributions in the photon. The need for a 
resolved photon interaction, i.e. where a photon interacts via its partonic agents, has become 
apparent in other type of processes involving photons, namely in the production of particles 
with a large transverse momentum. A recent review of the experimental situation and the 
existing parametrizations can be found in ^. 

Our motivation for a new global analysis of the data and for constructing a new parton 
parametrization for a real unpolarized photon is twofold. On the one hand, there is a vast 
amount of new experimental data on F2{x,Q'^) that has not been used yet to produce the 
)arton parametrizations for the photon. Two recent parametrizations, GRV and GRS 
^, used respectively about 70 and 130 experimental points, while at present a total of 208 
independent Q^) points exist. On the other hand, there are discrepancies between the 

theoretical calculations and experimental results for some processes initiated by real photons 



in which heavy quarks are produced|49|. Let us just mention here the D*- and D^-meson 
photoproduction 0, B and 3 or the Z}*-meson production with associated dijets [7[ at 
HERA, as examples. The disagreement between the theoretical and experimental results 
is even more profound for the open beauty-quark production in both HERA and LEP 
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12| | measurements. 



The idea of the radiatively generated parton distributions has been successfully intro- 
duced_by the GRV group first to describe the parton distributions in the nucleon [l3| and 
pion |il4|]. and later to create the LO and NLO parton parametrization for the real !j] and 
virtual jl5| photon. Here we follow this approach for a real photon case, limiting ourselves. 
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to the analysis based on the LO QCD. The NLO analysis is under preparation. 

As mentioned above there is a problem with the QCD description of heavy-quark pro- 
duction in processes initiated by photons. Therefore, our analysis especially focuses on the 
heavy-quark contributions to the F2{x,Q'^). We apply a new Variable Flavour-Number 
Scheme (VFNS) approach, denoted by ACOT(x), proposed for heavy-quark production in 
the ep collision ( "electroproduction" ) in Q|. For comparison we perform a standard Fixed 
Flavour-Number Scheme (FFNS) fit as well. Since these two approaches are based on very 
distinct schemes, and since they need different evolution programs, we will refer to them as 
to two models, CJKL (ACOT(x) type) and FFNSc jii-L models, respectively. 

Our paper is divided into six parts. In section 2 we describe various approaches including 
the ACOT(x) scheme applied to the production of heavy quarks in hadronic processes. 
Section 3 is devoted to the description of the F2 in LO QCD, paying special attention to 
an implementation of the ACOT(x) scheme in the calculation of the F2{x,Q'^). In section 
4 a description of the two global fits performed by us is given. In particular, we present the 
solutions of the DGLAP evolution in both models. We describe in detail the assumptions 
for the input parton densities. In the fifth section of the paper, the results of the global fits 
are discussed and compared with the data for the F2{x, Q^), and for the -^^(x, Q^) averaged 
over various x-regions. A comparison with LEP data for Fg^^ is presented in section 5 as 
well. The summary of the paper and an outlook of work in progress can be found in section 
6. Finally, in the appendix we give a simple parametrization of the CJKL (LO) parton 
distributions. 

II. VARIOUS SCHEMES FOR A DESCRIPTION OF HEAVY-QUARKS PRO- 
DUCTION: THE PROTON-TARGET CASE 

In this section we describe various schemes, which are used in the calculation of the 
heavy-quark production in hadronic processes. There exist two standard schemes. In the 
FFNS, the light quarks {u, d and s) and the heavy ones (c and h) are treated on a different 
footing. The light ones are treated as being massless, and together with the gluons, are the 
only partons in the proton. The massive charm- and beauty-quarks are produced in the 
hard subprocesses: they can only appear in the final state of the process. In the second 
scheme, the Zero- mass Variable Flavour- Number Scheme (ZVFNS), when the characteristic 
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hard scale of the process is larger than some threshold associated with a heavy quark, this 
quark is also considered as a massless parton in the proton, in addition to the three light 
quarks. In this way, the number of different types of quarks (flavours) that we treat as 
partons in the proton increases with the scale of the process. 

For the Deep Inelastic Scattering on the proton (DISep), where the structure functions 
of the proton are measured, the condition for considering a heavy quark to be a parton of 
the proton is given by a simple (kinematic) threshold condition for the total energy in the 
7*p collision W, namely W > 2mh (from now on we will denote the heavy quarks by h). It 
defines the kinematically allowed region for the production of a heavy-quark pair. However, 
for the structure functions, e.g. F^x, Q^), and further for the parton densities g^(x, Q^), 
not W but the virtuality of the probing photon is considered to be a natural scale. In 
the inclusive production of heavy quarks, their transverse momentum or mass is often taken 
as a characteristic (hard) scale /i, where /i ^ ^qcd, at which parton densities of the initial 
hadrons are probed. The two, massive and massless, approaches are considered to be reliable 
in different fi regions. The FFNS loses its descriptive power when fi ^ rrih', on the other 
hand the ZVFNS does not seem appropriate if /i ~ rrih. In order to achieve a prescription 
working in all hard scale regions, various schemes trying to combine the two approaches have 
been proposed. They have a generic name: Variable Flavour- Number Schemes (VFNSs). 
The first of such approaches was introduced by the ACOT group in [l7|. Oth er g roups, such 
as RT Q, BMSN or CSN Q, 

created their own versions of the VFNSl5flll. 

Let us discuss some aspects of the VFNS in more detail. The VFNS introduce the notion 
of "active quarks", for which the condition fi > niq is fulfilled. Such quarks can be treated 
as (massless) partons of the initial hadron(s). Light {u,d and s) quarks are always active 
because for them fi ^ Aqcd > "^g. The heavy-quark densities vanish for fi < rrih, otherwise 
they differ from zero like in the ZVFNS. For example, for the charm quark we see that at 
fi = riic we turn from a Three Flavour-Number Scheme to a Four Flavour-Number one. If 
Nf equals the number of active (massless) quarks, we define the Nf-FNS as one where the 
Nf first quarks are treated as light and the remaining quarks as heavy. 

In calculations based on the VFNS we take into account the sum of all contributions, 
which would be included separately in the ZVFNS and FFNS. Such procedure requires a 
proper subtraction of the double-counted contributions. Such contributions have the form 
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of the large logarithms In /i^, and are already resummed in the density qh{x, /i^). 

An important aspect of the VFNS is the behaviour of the heavy-quark contributions in 
the threshold region. Let us discuss as an example the production of heavy quarks in DlSgp. 
As was already mentioned, a heavy quark can be considered as a parton of the proton if 
the centre of mass energy of the hard process is W > 2mh. However, if we use /i^ equal 
to Q^, where = W^^^, and impose the threshold condition on /i^, then, for any Q^, it 
may happen (for small enough x) that qhi^iQ'^) = in the kinematically allowed region, 
i.e. for W > 2mh. On the other hand, non-zero heavy-quark densities may appear in the 
kinematically forbidden region in the (x, Q'^) plane. Moreover, such conditions can lead to a 
very steep or even non-continuous growth of the heavy-quark distributions at the threshold. 

In general one should ensure that all the ZVFNSs and relevant subtraction terms smoothly 
vanish when W —>■ 2mh. Then, the non-zero contributions should give only those terms 
which arise in the FFNS approach, since this approach should reliably describe the region 
W ~ 2mh. Different threshold conditions were used in different analyses; in particular, the 
ACOT group proposed to use a variable /x^ given by 

ml + cQ\l-ml/Q^Y for > ml 
= { (1) 
m\ for < ■'^L 



where c = 0.5 and n = 2 in 



17[. Still, in their approach the heavy-quark densities satisfy 



the boundary condition at Q = mj^ 

f=0 for ii^ = ml {Q'<ml), 
I 7^ for /i^ > ml (Q^ > ml). 

Recently, a purely kinematic solution of the threshold-behaviour problem has been found, 
on which the so-called ACOT(x) scheme is based. A new variable, Xh = xil + Aml/Q"^), 
has been introduced to replace the Bjorken x as an argument in the heavy-quark h density 
in the ZVFNS contributions. More details can be found in section 3. 

Although the above discussion was focused on the proton case, the problems with the 
proper treatment of the heavy-quark thresholds in the parton distribution are very similar 
for any target. In this paper we adopt the ACOT(x) scheme to the real photon case for the 
very first time. 
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III. DESCRIPTION OF IN THE ACOT(x) SCHEME 



In this section we first recall the basic facts related to the structure function -F^(x, Q^) 
for the real photon. Then we introduce the ACOT(x) approach for the photonic case. 



A. The parton densities in the photon 

The Deep Inelastic Scattering on a real photon (DlSg^) allows us to measure the structure 
function F^, and also other structure functions, F^, -F^, ■ ■ ■ , via the process 

7*7 — i> hadrons, (3) 

see [1-3]. In the Parton Model this is described at lowest order by the Bethe-Heitler (BH) 
process, 7*7 — > qq (see Fig. [T]). 

In the leading logarithmic (InQ^) approximation or, in short, in the leading order of 
QCD (LO QCD), the photon structure function Q"^) can be written in terms of quark 

(antiquark) densities qj{qj) as follows 

1 

-f:1{x,Q') = Y.e'Ml + (4) 

X 

1=1 

where Nf is the number of different quark flavours, that can appear in the photon ( "active 
quarks"). Note that qJ{x,Q'^) = q'J{x,Q^). 

The evolution of the parton densities with InQ^ is governed by the inhomogeneous 
DGLAP equations. In LO we have for a quark (similarly for an antiquark) and a gluon 
density 



dq]{x,Q^) a 2, / ^ , ^siQ^) dy 
— —e^k[xj + ' 



dlnQ2 27r 



Nf 

i^A E(gZ + ?7)(l/, Q') + Pgg i^A G\y, Q') 



(5) 
• (6) 



The k{x) term on the right-hand side of Eq. ^ comes from the Bethe-Heitler process 
of Fig. m for 3 colours we have 

k{x) = 3 [x^ + {1 - xf] . (7) 
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FIG. 1: The Bethe-Heitler process 7*7 — > qq. 



The functions Pi{x) are the LO sphtting functions 22 
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^.G(x) = -[a;2 + (l-xf] 
- - ^ , 



Pgg{x) = 6 



H — + a;(l — x) 

X 





'11 


Nf 


+ 




_Y 


3 



(8) 



5(1 - a;). 



L(l-x)+ 

Note that the function k{x) describes a photon into quark sphtting, so one has k{x) 



Pg^{x). 



1. Heavy-quark contributions to in the FFNS 

The standard FFNS approach corresponds to a number of "active quarks" (Nf) =3, so 
only the hght quarks (and their antiquarks) are taken into account in the sum in Eq. (^. The 
main heavy-quark contributions to are obtained in this scheme from the corresponding 



Bethe-Heitler process (Fig. 1 with q —>■ h) 

7*7 hh, 

keeping the heavy-quark masses in the calculation. It reads 

1 ry 
-Fl,ix,Q')U, = 3-eiu;{x,Q') 



(9) 



(10) 
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with 



1 + 8x(l — x) — xil — x) 



Ami 



+ ln 



1 + P 
1-/3 



x^ + {l-xf + x{l-3x)^ -x^ ^ 



Amlx 



1 - 



Ami 



(11) 
(12) 



(l-a;)g2 V W^' 

We call this contribution -F2^^|dir since here the real photon (i.e. the target photon) interacts 
directly. However, there exists another heavy-quark contribution, related this time to the 
process with the resolved initial photon, namely 

7*G^ ^ hh, (13) 

with a gluonic parton of the photon target (as in Fig. Qwith 7 G^); it gives 

F^^,{x,Q')U = ^^el y ^c.(^,Q2) G^{z,Q')dz, (14) 



'Xh 



where 



Xh = x{l + 



Ami 



(15) 



2. Heavy-quark contributions to F2 in the ZVFNS 

In the ZVFNS, the number of "active quarks" changes with the hard scale, as described 
in the previous section. For low scales the sum in Eq. extends to Nf = 3 but whenever a 
heavy quark threshold is surpassed the value of Nf is increased by 1. It is worth mentioning 
that in some parton parametrizations for a real photon the heavy-quark densities do appear; 
however they are described in the threshold region W ~ 2mh only by the above Bethe-Heitler 
formula (Eqs. (10) and (11)). Moreover, instead of the restriction on W one sometimes 
takes a (reasonable) condition on (see, for instance, the GS parametrization !3]). Of 
course well above a heavy-quark threshold, such a quark can be included among the active 
(massless) quarks and then Nf — > Nf + 1, see e.g. 14\- 



B. ACOT(x) scheme for 



The ACOT(x) prescription combines the FFNS and ZVFNS, so that we have to add all 
relevant contributions from both approaches. For the light-quark contributions we take the 
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form given in Eq. Q with Nf = 3, while for the heavy quarks we include the following 
terms: 

F^{x,Q')U,= [xe^(g,^ + g;:)(a:,Q2)+F,'^,(x,Q^)U. + F2:.(^,Q')Us], (16) 

h{=c,b) 

where F2''^(x, (5^)|dir and (5^)|res are given in Eqs. (fTUI) and (fT^ . respectively. 

In Eq. (fTBj) we double-count some heavy-quark contributions. Indeed, part of the F^f^Air 
contribution from 7*7 — > hh corresponds to the coUinear configuration. Such a configuration 
leads to a contribution proportional to In and is already included in the DGLAP equation 
for ql{x, Q^), via the k{x) term. Therefore we must subtract from (|T6|) the following terms 

F2^hUn,sMhtv = X In -^3e^- (x^ + (1 - xf) , (17) 

IT 

coming from an exact solution of a part of the DGLAP equation, namely 



dql{x,Q^) _ a 



eiHx), (18) 



rflnQ2 2tt 

integrated over the from ml to Q^. 

Similarly, Fg'^^lres from the 7*G'^ — hh process, has a InQ^ part that corresponds to 
the coUinear configuration already included in the DGLAP equation for gj[(a;, Q^), via the 
PqG j G^{y, Q^) term. The term to be subtracted reads, in this case: 



F2 h I reSjSubtr X 



\^hJ ^ Jx y \yj 



It is based on an approximated solution for the other part of Eq. (jSl), namely 
dql{x,Q^) _ a,(Q2) fUy 



X \ , 



(20) 



dlnQ^ 27r J.^ y 

The solution (|Tn|) is obtained by the inte gra tion of Eq. ()2()|1 over the same region as 
above, after neglecting the dependencejMi of as and G^. The final subtraction term is 
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F2 {X-iQ )|subtr;c,6 — ^ ^ [-^2,h I dir,subtr ~l~ -f2,/i I res,subtr] • (^l) 
h{=c,b) 

So, finally we have, for heavy-quarks: F^{x,Q'^)\c,h = F^{x,Q'^)\c,b - F^{x,Q'^)\suhtr,c,b- A 
graphical representation of all terms included in the analysis, Eqs. (fT^ and (j2H), is presented 
in Fig. El 




FIG. 2: Graphical representation of the ACOT model for F2 ■ The first diagram represents the 
ZVFNS contribution. The third (fifth) diagram shows the FFNS contribution of the resolved 
(direct) photon, while the second (fourth) diagram is the corresponding subtraction term. 

Further, we need to ensure that all terms containing the heavy-quark h disappear when 
W 2mh. The FFNS contributions, F2^;,|dir (Eqs. (UMl) and F2^;,|res (Eqs. P|).(fT3|l). 
behave properly at the thresholds. The problem emerges for the heavy-quark densities 
and the subtraction terms (Eqs. (fTH|) . (fT^ ). These terms do not naturally disappear for 
W — i> 2m/j. Fortunately, this problem can be cured; we noticed that the resolved-photon 
contribution in Eq. ()14|) vanishes for W —>■ 2mh because then Xh —>■ ^ and the corresponding 
integral disappears. So, we can do the same with the q1 distribution and the subtraction 
terms, if instead of x we introduce the Xh variable (Eq. (fTH|l ) slightly shifted from x. 
This way we force the heavy-quark distribution and the second term of the subtraction 
contribution (the integral term) to vanish at the corresponding threshold. Unfortunately, 
unlike for the proton, in the case of the photon we are left with the -F2^/j|dir,subtr contribution, 
which is now proportional to + (1 — Xh)'^ and does not vanish for Xh 1- In the 
large-Q^ region, where the ZVFNS is reliable, this change of variables is irrelevant. In the 
numerical calculations we ensure that the total contributions to F2 due to heavy quarks 
are not negative (positivity constraint). This way we effectively introduce small additional 
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terms near the charm- and beauty-quark thresholds. The final formula for the F2{x, Q"^) in 
the ACOT(x) scheme reads {52! 

F^x^Q') = Y.xei{qJ + q]){x,Q')+ ^^Hql + qDixn.Q') 

i=l h{=c,b) 
2 

+ E [FU^,Q')Ur + Fl,{x,Q')U] (22) 

h{=c,b) 

h{=c,b) -^Xh^ \y/ 

As can be seen the heavy-quark distributions are included in the second sum as qKxh, Q'^), 
being Xh functions of x (and Q^). We parametrize the final form obtained for these distri- 
butions in the Appendix as simple functions of x and Q^. js^ 

The Xh variables of the ACOT(x) scheme recall the so-called "slow rescaling" obtained 
in early papers on the charm-quark production in the DISep (e.g. j2^), where the Bjorken 
X was replaced with (: 

IV. GLOBAL FITS - SOLVING THE DGLAP EVOLUTION 

Using all the existing F2{x,Q'^) data we perform two global fits, both based on the LO 
DGLAP evolution equations. One fit uses the ACOT(x), the other a FFNS model for the 
heavy-quark contributions. 

First we introduce the basic ingredients that are common for the two considered models. 

A. Mellin moments 

The LO DGLAP evolution equations are very much simplified if they are transformed 
into the Mellin-moments space. The n-th moment for the quark or gluon densities, is 
defined by 

' (24) 



Jo 
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Analogous definitions can be used for tlie splitting functions PJj. The evolution equations 
in the Mellin space have the generic form 



(25) 



dlnQ2 

Obviously, the first term on the right-hand side appears only for the quark densities. For 
simplicity, in the following we will skip all subscripts and superscripts wherever possible. 



B. Non-singlet and singlet parton densities 

To solve the DGLAP equations we need to decompose the parton densities into the singlet 
and non-singlet (in flavour space) combinations. For the non-singlet (ns) case we have 



i=l 



(26) 



^ = 2Ar^ ( (e^) - {e'f) k 



where 



Nf 



(27) 



i=l 



and Cj stands for the corresponding quark electric charge. Similarly for the singlet (s) 
densities we have 



s^(g2) 



(28) 



i=l 

In the singlet case the DGLAP equations ()25p become a matrix equation with 



p = 


p 

J- qq 


^NfPqG 


k = 


ks 




Pcq 


Pgg 








with 



h = 2Nf{e^)k. 



(29) 



(30) 
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C. Point- and hadron-like parts 



The solution of the DGLAP equations can be divided into the so-called point-like (pi) 
part, related to a special solution of the full inhomogeneous equation and hadron-like (had) 
part, arising as a general solution of the homogeneous equation. Their sum gives the partonic 
density in the photon, so we have 



rm = fUQ') + fW) (31) 



where 



pfn"^) - 1 " ri _ rl-2P//3ol /, 

(32) 

/had(g^) = L-^^/*r(g^). 

Here jSo = 11 — 2Nf/3, P equals Pgg (P) for the non-singlet (singlet) parton densities, and 
L = "^Ig^jj , where the Qq is the evolution starting (input) scale. Note that at = Qq the 
point-like part vanishes. 



D. Input parton densities. VMD 

Following 0|, the input scale has been chosen to be small, Qq = 0.25 GeV^, hence our 
parton densities are radiatively generated. The point-like contributions are given by Eq. 
(30), while the hadronic parts need the input distributions. For this purpose we utilize the 
Vector Meson Dominance (VMD) model j25|, where 

/h\d(^,Q^) = E^/''(^'^o), (33) 
V Jv 

with the sum running over all light vector mesons (V) into which the photon can fluctuate. 
The parameters fy can be extracted from the experimental data on the T{y — > e^e~) width. 
In practice one takes into account the p*^ meson 5J] and the contributions from other mesons 
are accounted for via a parameter which is left as a free parameter. We take 

fLMl) = ^^rMl)- (34) 
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In the GRV prescription 



a 



the parton densities in the meson are approximated by 
the pionic ones: /^(a;, Ql) ~ /'^(x, Qq)- However, this assumption ignores among others the 
possible effects of the pseudo-Goldstone boson nature of the pion, and we are not using it in 
our analysisQl- Instead we use the input densities of the p° meson at = 0.25 GeV^ in the 
form of valence-like distributions both for the (light) quark (v^) and gluon {G'') densities. 
All sea-quark distributions (denoted by are neglected at the input scale. At this scale, 
the densities v^, and are related, according to Eq. (jMj) to the corresponding densities 
for a photon; see below. 
The density is given by 

v'{x, Ql) = ^{u^^ + uP~ + dP~ + dP^){x, Ql), (35) 

where from the isospin symmetry 

u"^ {x,Ql) = u'' (x, Ql) = (x, Ql) = dP^ {x,Ql). (36) 

Note that all the densities in Eq. (35) are normalized to 1, i.e. u^^dx = 1. 

The following constraints should hold for the f ''(x, Ql) density. The first is related to a 
number of valence quarks in the p° meson, and we have 

»i 

2vP{x,Ql)dx = 2. (37) 



The second constraint represents the energy-momentum sum rule 



/ x{2vP{x,Ql) + GP{x,Ql))dx=l. (38) 
Jo 

We parametrize the input densities as follows 

xCP{x,Ql) = 0, (39) 
xvP{x,Ql) = NyX°{l-x)^, (40) 
xGP{x, Ql) = Ngxv^ix, Ql) = Ngx'^il - x)^, (41) 

where Ng = NgNy, and impose two constraints given by Eqs. ()37|) and (jHHj) in both models. 
These constraints allow us to express the normalization factors A^^, and Ng as functions of 
a, (3 and k. This leaves these three parameters as the only free parameters to be fixed in 
the fits to the F2 experimental data. 
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E. FFlSiScjKL model 



In the FFNScjKL model the number of "active quarks" (Nf) equals 3, so (e^) = 2/9. 

To describe the hadron-like part of the solution of DGLAP equations for the photon, we 
introduce, apart from the valence-like quark and gluon densities, also the sea distribution 
Chad(^; Q^)- "^^^ ^P"; down- and strange-quark densities in the photon are then given by the 
following combinations 

<d(^, Q') = dUi^, Q') = \ Kd(^> Q') + 2Cad(^, Q')] , (42) 
sUM') = CU^,Q'). (43) 
/^From Eqs. (26) and we get (below we simply use f2s instead of /Jsg): 

/nlhad(^,Q') =l<A^,Q'), (44) 

sLd(^, Q') = Q') + 6Cad(^> Q')- (45) 

After performing the DGLAP evolution of /nshad(^'Q^) ^had(^'Q^) fro^^ Qo to higher 
Q^, we calculate v^^^{x, Q"^) and Chad(^' Q^)- Finally, using formulae P2|l and (j43|l . we obtain 
the hadron-like part for the individual quark densities. 

As the down- and strange-quarks have equal electric charges, there are only two different 
point-like distributions: Upi{x,Q^) and dp[(a;, Q^) = Spi(a;,Q^). We calculate them again 
through the evolution of the singlet and non-singlet combinations of the parton densities. It 
can be easily checked that distributions read as 

1 (46) 

Finally, the contribution due to the massive c- and 6-quarks are approximated by the 
Bethe-Heitler formula (10)-(12) for Q"^ > Ami. 

F. CJKL model 

In the CJKL model all terms originating from both FFNS and ZVFNS are included. This 
means that apart from the light-quark distributions we take into consideration c'^{x, Q"^) and 

15 



W{x, Q^), which emerge in the ZVFNS, so here Nf = 5. 

When five "active quarks" are considered instead of three, the DGLAP evolution becomes 
shghtly more comphcated and we need more non-singlet parton densities than for the simple 
FFNS model. Here we need /nsa, /nss, /ns4, and /nsg (Eq. (UHl)) for both hadron- and point-like 
parts; using them and S, calculated for Nf = 5, we can express individual quark densities 3] 



as follows: 



1 



P = ^i+mi, - 45/4 - 30/n\ + 15/4 + 2S^), (47) 
c'' = ^( - 45/4 + 30/4 + 15/4 + 2S^), 

b' = ^i +30/4-30/4+ S^). 

For the hadron-like parts we consider, similarly to the FFNS case, the light-quarks den- 
sities given by Eq. (021) • Among sea quarks we have now all types of quarks, in particular 
we have: 

s'U^, Q') = c'Ux, Q') = b]ax, Q') = CU^, Q'). (48) 
This leads to the following relations 

p = p = n 

Jns2,had Jns4,had 

/ns3,had ~ "^"^had' (49) 

P = —iP 

•'nssjhad 15 

valid at every x and Q^. 

In the point-like case equality of the electric charges for the up-type and down-type 
quarks leads to the following relations: u^^{x, Q"^) = Cpj(x, Q"^) and d^i{x, Q^) = s'^^{x, Q^) = 

In our analysis, performed within the ACOT(x) scheme, we evolve first the singlet and 
non-singlet distributions and we obtain in this way the gluon g{x,Q'^) and light-quarks 
M, d, s(x, Q"^) densities. In each heavy-quark density, i.e. for c and h quarks, we replace x 
variable by the corresponding Xh variable; let us recall that Xh = x\\^ \ , with mh equal 
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to rric and rrib, respectively. For such densities we perform DGLAP evolutions and obtain 
the resulting c{xh,Q^) and b{xh,Q^) distributions. We then compute the F2{x,Q^) using 
Eq. and fit the parameters of the model to the experimental data, which is described 
in more detail in the next section. 

Finally, we parametrize all our resulting parton distributions analytically with the possi- 
bly simple functions of x and Q^. This parametrization is given in the Appendix. 



V. GLOBAL FITS AND RESULTS 



All together, 208 



We have performed two fits to all the available Q ) data j28l|- 

data points were used, including the recent high-Q^ measurement of the OPAL collaboration 
Isoj ]. Some of the data points are not in agreement with others. We will discuss in detail 
their influence on the fit in the next section. The fits based on the least-squares principle 
(minimum of x^) were done using MiNUlT Systematic and statistical errors on data 
points were added in quadrature. 

The as value used in the fits was calculated from the LO formula, which depends on Nf 



For Nf = Awe took the QCD scale A^^) equal to 280 MeV j41] with the assumption that the 

LO and NLO A values for four active flavours are equal, which is consistent with the GRV 
I I 

group approach [4q. Values of A for other Nf can be calculated if one assumes a continuity of 
as at the heavy-quark thresholds = m|. Assuming then that as{mfi)^^f^ = as(m|)^^^+^^ 
one obtains the relation 

A(^/+i) = m,(A(^/)/m,)(=^'-'^^)/^=^'-'^^\ (51) 

which gives A(3) = 314 MeV and A^^) = 221 MeV. Finally, the masses of the heavy quarks 
are taken to be 41]: rric = 1.3 GeV and mj, = 4.3 GeV. 

The results of both fits are presented in Table HI The second and third columns show the 
quality of the fits, i.e. the total for 208 points and the per degree of freedom. The 
fitted values for parameters a, (3 and k are presented in the middle of the table. In addition, 
the values for N^ and Ng obtained from these parameters using the constraints (|n7j) and 
are given in the last two columns. 

17 



Model 




X^/dof 


K 


a 


/3 






FFNScj/^L 


All 


2.30 


1.726 


0.465 


0.127 


0.504 


1.384 


CJKL 


430 


2.10 


1.097 


0.876 


2.403 


2.644 


2.882 



TABLE I: The for 208 points and parameters of the fits for Y'F'HScjkl and CJKL models. 



We see that the obtained per degree of freedom is better in the CJKL model than 
in the standard-type FFNS approach; however, it is not particularly good, owing the poor 
quahty of some data used in the analyses. This fact has been already discussed in many 
papers, e.g. |4jJ], see also discussion in section 5.1. 

The two fits to the same collection of data, although not very different as far as is 
concerned, are obtained with very different sets of parameters. Note that k is close to 1 in 
the CJKL case, while for the other fit it is closer to 2. If this parameter is close to 1, we have 
in practice only the p contribution at the input scale. However this is not the whole story 
since the nN^ and nNg give full normalizations of the valence-like quark and gluon densities 
in the 7. Now, the and A^^ are much smaller in the standard approach than in the CJKL 
model. Finally, let us notice the large difference in both models, small and large x, of the 
fitted input densities, which correspond to very different a and (3 parameters, respectively. 
The FFNSc model has a close to the standard (Regge) one for a valence-quark density 
(a — 1 ~ —0.5); however its /5, which governs the large-x behaviour, is very small, being 
far from 2, a standard prediction from the quark- counting rule j^^]. On the other hand, the 
CJKL model gives j3 closer to 2, while its a — 1 is close to zero. 



A. Comparison of the CJKL and YY^^cjkl fits with F2 data and other LO 
parametrizations 

The x^ values, representing the quality of our LO fits, are compared in Table HTl with the 
corresponding x^ obtained by us using the GRS LO 

describe the present F2 data. The comparison is performed for a set of 205 data points, 
i.e. excluding the points with < 0.26 GeV^ since they were not used in performing the 
GRS parametrization. The second column gives the number of independent parameters in 
each model. The overall x^ and x^/DOF values are given in the middle of the table for 
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205 data points. It is clear that both our fits give a better description of the experimental 
data than the previous parametrizations. This could be expected since we are including 
more data in our fits. The CJKL model gives the lowest value of x^/DOF, but it is still 
rather high. This may arise from the fact that we use all available data and, as it was stated 
(e.g. j43|]), the data published by the TPC/27 Collaboration [3]j are inconsistent with other 
measurements. We studied this point and in the last two columns of the table we present 
the values calculated without the TPC/27 points. Indeed the x^/DOF then computed 
is visibly improved [s?!. A special CJKL fit performed without TPC/27 data gives x^/DOF 
equal to 1.78. The very recent NLO analysis performed in [4^ for 134 experimental points 
and with five independent parameters gives x^/DOF = 0.93. 



model 


# of ind. par. 


# of data points 


205 


182 - no TPC 




xVdof 




xVdof 


SaSlD 


6 


657 


3.30 


611 


3.47 


GRS LO 





499 


2.43 


366 


2.01 


FFNScji^L 


3 


442 


2.19 


357 


1.99 


CJKL 


3 


406 


2.01 


323 


1.80 



TABLE II: Comparison of values obtained for the FENScje-l and CJKL fits to the F^ix^Q'^) 
data with calculated using the SaSlD and GRS parametrizations. 

Figures OHEl show a comparison of the CJKL and FFNScji^L fits to F2{x,Q^) with the 
experimental data as a function of x, for different values of Q^. Also a comparison with 
the GRS LO and SaSlD parametrizations is shown. (If a few values of are displayed 
in a panel, the average of the smallest and biggest one was taken in the computation.) As 
can be seen in Figs. ElandlU both CJKL and FFNScjft-L models predict a much steeper 
behaviour of the -F7(a;, Q'^) at small x with respect to other parametrizations. In the region 
of a; > 0.1, the behaviour of the -F^(x, Q^) obtained from the FFNScji^/, fit is similar to the 
ones predicted by the GRS LO and SaSlD parametrizations. The CJKL model gives lower 
prediction whenever the charm-quark threshold is surpassed, and slightly below this region, 
as is shown in Figs. ElandlHl 
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Apart from this direct comparison with the photon structure-function data, we perform 
another comparison, this time with LEP data that were not used directly in our analysis. 
Figures □ and e present the predietions for F^i..Q% averaged over var-ous . regions, 
compared with the recent OPAL data 3!/|. For comparison, the results from the GRS LO 
and SaSlD parametrizations are shown as well. We observe that all models give very similar 
predictions, which are in fairly good agreement with the experimental data. Only the CJKL 
model slightly differs from the other models considered: it gives better agreement with the 
data. The difference between the predictions of the CJKL model and the other models is 
most striking in the case of the medium-x range, 0.1 < x < 0.6 shown in Fig. The CJKL 
curve clearly shows a departure from the simple InQ^ dependence. This is caused by the 
additional dependence due to the x variable. The highest x range, 0.85 < x < 0.98 (see 
Fig. IHl) is the second region of a significant difference between the predictions of the CJKL 
and other models. The predictions split at high-Q^ values, as the CJKL predicts a much 
softer dependence. 



B. Parton densities 



In this section we present the parton densities obtained from the CJKL and FFNScje-l 
fits and compare them with the corresponding distributions of the GRV LO ^, GRS LO 
and SaSlD parametrizations. First, we present all parton densities at = 10 GeV^ (Fig- 
E}. The biggest difference between our CJKL model and others is observed, as expected, 
for the heavy-quark distributions. Unlike for the GRV LO and SaSlD parametrizations, the 
densities c'*'(x, Q^) and W{x, Q"^) vanish not at x = 1 but, as it should be, at the kinematical 
threshold. Also the up-quark density differs among models. In the CJKL model it is lower 
than in other parametrizations for x > 0.1. The same can be seen in Fig. ^1 where 
for various values the up-quark distributions are presented. The up-quark density in the 
CJKL, FFNScjKL and GRV LO models have similar behaviour at very small x. The hardest 
up-quark distribution is obtained in the FFNSc ji^^ approach, while both the GRS LO and 
SaSlD predictions are much softer. The same holds in the case of the gluon distribution, 
shown in Fig. ^2 Finally, the charm-quark densities of the CJKL model and of the GRV 
LO and GRS LO parametrizations are presented in Fig. ^| Here we see that, in addition 
to the already mentioned vanishing at the threshold, the charm-quark distribution in the 
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CJKL model is larger than the ones in the other parametrizations. This is particularly true 
for larger values of Q^, where the threshold is very close to a; = 1. 

Finally, in Fig. we present our predictions for the F^^. For = 5, 20, 100 and 
1000 GeV^ we compare the individual contributions included in the CJKL model. As was 
explained in detail in section 3.2, they all, apart from the -F2,c|subtr,dir term, vanish in the 
W — i> 2mc threshold. This term dominates near the highest kinematically allowed x. The 
direct BH term is important in the medium- a; range. Its shape resembles the valence- 
type distribution. The charm-quark density contribution, i.e. the term 2xe'^c'{x,Q'^), is 
important in the whole kinematically available x range; it dominates the for small x. In 
this region also the resolved-photon contributions increase, but they cancel each other. 

C. Comparison with 

A good test of the charm-quark contributions is provided by the OPAL measurement of 
;he obtained from the inclusive production of D*^ mesons in photon-photon collisions 
4^ . The averaged has been determined in the two x bins. These data points are 
compared to the predictions of the CJKL and FFNSc ji^L models as well as with the SaSlD 
and GRS LO parametrizations in Fig. The prediction of the FFNScji<:L niodel, which 
as the only one among those compared does not contain the resolved-photon contribution, 
is based on the point-like (Bethe-Heitler) contribution for heavy quarks only. As seen in the 
figure it decreases too quickly with the decreasing x, much faster than the predictions from 
other parametrizations. The CJKL model seems to give the best description of the data for 
the low-x bin, but overshoots the experimental point at high x. 

VI. SUMMARY AND OUTLOOK 

A new analysis of the radiatively generated parton distributions in the real photon based 
on the LO DGLAP equations is presented. All available experimental data have been used 
to perform two global, 3-parameter fits. Our main model (CJKL) is based on a new variable 
flavour-number scheme (ACOT(x)), applied to the photon case for the very first time. 
It has a proper threshold behaviour of the heavy-quark contributions. The CJKL-model 
results are compared with an updated Fixed Flavour-Number Scheme (FFNScji<:L) fit and 
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to predictions of the GRS LO and SaSlD parametrizations. Our model gives the best 
of those compared. It describes very well the evolution of the F^ix^ Q"^) , averaged over 
various x-regions. We have checked that the CJKL fi t ag rees also reasonably well with the 
prediction of a sum rule for the photon, described in 45 1. 

We have also checked that the gluon densities of both CJKL and FFNSc jt^l niodels 
agree with the HI measurement of the gluon density (C^) performed at = 74 GeV^ 
|4q . In both models gluon densities are very similar to the gluon density provided by the 
GRV LO parametrization, which gave so far the best agreement with the HI data. Further 
comparison of our gluon densities to the HI data cannot be performed in a fully consistent 



way 



58l | , since GRV LO proton and photon parametrization were used in the experiment in 



order to extract such gluon density. 

One of the motivations for this work was given by the disagreement between the theo- 
retical and experimental results for the open beauty-quark production in two-photon pro- 



cesses in the LEP 



12| measurements. We did calculate the LO cross-section for charm- 
and beauty-quark production in 77 collisions in the ACOT(x) and FFNS schemes, us- 
ing the CJKL and FFNSc jkl distributions of partons, respectively. The cross-section for 
the c-quark production computed in both models agrees with the experimental data. The 
ACOT(x) model gives a slightly better shape of the c-quark distribution. There is a small 
difference between the results of the two models for 6-quark production. We observe an in- 
crease of the cross-section for the beauty-quark in the ACOT(x) approach, as compared to 
the FFNS result, based on GRS LO parametrization, but it is too small to fit the experimen- 
tal data. More work on this subject is required. Also, before reaching a definitive conclusion 
the NLO corrections should be considered and the NLO parton densities for the photon ap- 
plied. The NLO parametrization in the CJKL model, together with the effects of different 
subtraction terms and positivity constraints, will be presented in a future publication. 

A simple analytic parametrization of the results of our CJKL model for the individual 
parton densities is given, and a FORTRAN program calculating parton densities as well as 
the structure function F2 (Eq. (jSH))) can be obtained from the web sites given in 
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APPENDIX: PARTON PARAMETRIZATION FOR THE CJKL MODEL 

We give here an analytic form for the parametrization of the CJKL results for the indi- 
vidual parton densities. Following the GRV group we parametrize them in terms of 

ln[Q7(0.221 GeV)^] 

ln[Q2/(o.221 GeV)2] ^ 

with Qq = 0.25 GeV^. The parametrization has been performed in the 10~^ < a; < 1 
and 1 < < 2 X 10^ GeV^ ranges. We made a separate parametrization of the point- 
and hadron-like densities. The parametrized distributions of the light {u, d, s) quarks are 
in agreement with the ones obtained in the fit up to few percent accuracy. The same 
is true for the gluon density, apart from the high-x region (where values fall down), 
where 10% accuracy is assured. We checked that the heavy-quark densities, for c- and 
6-quarks, are represented by our parametrization at 10% accuracy for > 2 GeV^ and 
10 GcV^ respectively. For both quarks those vahics are below their masses squared 
= 1.3^ and 4.3^ GeV^, which are often considered as the energy scale of the processes 
involving heavy quarks. The obtained for 205 data points for Q^) is equal 406, for 

both fitted and parametrized distributions. 
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The final parton densities in the real photon are (we skip the superscript 7 in this part): 

G(x,g2) = Gpi(a;,Q2) + Ghad(a;,Q') 

d{x,Q^) = dp,{x,Q'') + ^vix,Q')+C{x,Q'') 

s{x,Q^) = Spi{x,Q^)+C{x,Q^) (A.2) 

c{x,Q^) = Cpi(a;,(5^) + Chad(x,(5^) 

b{x,Q^) = 6pi(x,g2) + 6had(a;,Q^) 

Note that all our densities describe the massless partons, although the kinematical con- 
straints for c- and 6-quarks were taken into account. The formulae given in the Appendix 
parametrizing densities of both heavy quarks, represent the densities as included in the sec- 
ond sum in the Eq. (j2Sl), that means that, for instance, the final density b{x, Q"^) should be 
understood as being equivalent to b{xh{x,Q'^),Q'^). 

A FORTRAN code of the parametrization can be obtained from the web pages j^]. The 
program includes also an option for the F2{x, Q"^) calculation according to the Eq. (|^. 

1. Point-like part 

We use the following form for the parametrization of the point-like distribution for light 
quarks and gluons (denoted by /pi(a;,(5^)) 



^"^^^^"'^'^ = 1^^" (0.22fGeV)^ 



X 



(A.3) 



s"x\A + + Cx^) + s"' exp {^E + \jE'sP In ^ j 
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For the gluon density Gpi: 



a = -0.43865 , a' = 2.7174 , /3 = 0.36752, 
A = 0.086893 - 0.34992s, 

B = 0.010556 + 0.049525s, , , 

(A.4) 

C = -0.099005 + 0.34830s, D = 1.0648 + 0.143421s, 

E = 3.6717 + 2.5071s, E' = 2.1944 + 1.9358s, 

a = 0.23679 - 0.11849s, b = -0.19994 + 0.028124s. 

For the up-quark density Upi: 

a = -1.0711 , a' = 3.1320 , P = 0.69243, 
A = -0.058266 + 0.20506s, 
B = 0.0097377 - 0.10617s, 

C = -0.0068345 + 0.15211s, (A.5) 

D = 0.22297 + 0.013567s, 

E = 6.4289 + 2.2802s, E' = 1.7302 + 0.76997s, 

a = 0.87940 - 0.110241s, b = 2.6878 - 0.040252s. 

For the down- and strange-quark densities, dpi — Spi: 

a = -1.1357 , a' = 3.1187 , P = 0.66290, 
A = 0.098814 - 0.067300s, 
B = -0.092892 + 0.049949s, 

C = -0.0066140 + 0.020427s, (A.6) 

D = -0.31385 - 0.0037558s, 

E = 6.4671 + 2.2834s, E' = 1.6996 + 0.84262s, 

a = 11.777 + 0.034760s, b = -11.124 - 0.20135s. 

In the case of heavy quarks a shghtly modified function hpi was apphed: 
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2 



l-KM')-l^- ^,^^^^^^y - (A.7) 



+ + Cy'') + s"' exp {-E+ ^^E'snn^ 



with y ^ x + 1 - Q2+6?6GeV^ charm-quark and y = x + 1 - Q2+73^96Gev^ 

beauty-quark densities. 

For the charm-quark density Cpi, for < 10 GeV^: 

a = 2.9808 , a' = 28.682 , /3 = 2.4863, 

A = -0.18826 + 0.13565s, B = 0.18508 - 0.11764s, 

C = -0.0014153 - 0.011510s, 

D = -0.48961 + 0.18810s, (A.8) 

E = 0.20911 - 2.8544s + 14.256s^ 

E' = 2.7644 + 0.93717s, a = -7.6307 + 5.6807s, 

b = 394.58 - 541.82s + 200.82s^ 

For the charm-quark density Cpi, for > 10 GeV^: 

a = -1.8095 , a' = 7.9399 , /3 = 0.041563, 

A = -0.54831 + 0.33412s, B = 0.19484 + 0.041562s, 

C = -0.39046 + 0.37194s, D = 0.12717 + 0.059280s, (A.9) 
E = 8.7191 + 3.0194s, E' = 4.2616 + 0.73993s, 
a = -0.30307 + 0.29430s, b = 7.2383 - 1.5995s. 
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For the beauty-quark density &pi, for < 100 GeV^: 

a = 2.2849 , a' = 6.0408 , (5 = -0.11577, 

A = -0.26971 + 0.17942s, 

B = 0.27033 - 0.18358s + 0.0061059s2, 

C = 0.0022862 - 0.0016837s, (A.IO) 

D = 0.30807 - 0.10490s, E = 14.812 - 1.2977s, 

E' = 1.7148 + 2.3532s + 0.053734^5, 

a = 3.8140 - 1.0514s, h = 2.2292 + 20.194s 

For the beauty-quark density 6pi, for > 100 GeV^: 

a = -5.0607 , a' = 16.590 , (5 = 0.87190, 

A = -0.72790 + 0.36549s, B = -0.62903 + 0.56817s, 

C - -2.4467 + 1.6783s, D = 0.56575 - 0.19120s, (A.ll) 

E = 1.4687 + 9.6071s, E' = 1.1706 + 0.99674s, 

a = -0.084651 - 0.083206s, b = 9.6036 - 3.4864s. 

2. Hadron-like part 

We use a simple formula for the valence-quark density: 

-xv(x,Q^) ^Ax''(l + By/^+Cx)(l-x)^, (A.12) 

a 

with the following parameters: 

A = 1.0898 + 0.38087s, B = 0.42654 - 1.2128s, 

C = -1.6576 + 1.7075s, D = 0.96155 + 1.8441s, (A.13) 
a = 0.78391 - 0.068720s. 



For the gluon distribution, we apply 



27 



a 

with 

a = 0.59945 , p = 1.1285, 
A = -0.19898 + 0.57414s, B = 1.9942 - 1.8306s, 
C = -1.9848 + 1.4136s, D = 0.21294 + 2.7450s, 
E = 1.2287 + 2.4447s, = 4.9230 + 0.18526s, 
a = -0.34948 + 0.47058s, b = 1.0012 + 0.99767s. 



In the case of the sea-quark density, we use: 

^xC{x,Q^) = + A^/^ + Bx) X 

X 

exp (^-E + ^E'sHn^ (1 - x)^, 

with 

a = 0.71660 , p = 1.0497, 
A = 0.60478 + 0.036160s, B = 4.2106 - 0.85835s, 
D = 4.1494 + 0.34866s, E = 4.5179 + 1.9219s, 
E' = 5.2812 - 0.15200s, a = 0.72289 - 0.21562s. 



Finally, for the heavy-quark densities: 



1 s° 

-xhh^{x,Q^) = + Ay/^ + By) X 



exp (^-E + )jE'snn^ (1 - y)^, 



28 



with y^x + 1- Q2+6^6GeV^ charm- and y^x + 1- Q2+r3.96Gey^ beauty-quark. 
For the charm-quark density Chad, for < 10 GeV^: 

a = 5.6729 , /3 = 1.4575, 

A = -2586.4 + 1910.1s, B = 2695.0 - 1688.2s, , . , 

(A.19) 

D = 1.5146 + 3.1028s, E = -3.9185 + 11.738s, 
E' = 3.6126 - 1.0291s, a = 1.6248 - 0.70433s. 

For the charm-quark density Chad, for > 10 GeV^: 

a = -1.6470 , P = 0.72738, 

A = -2.0561 + 0.75576s, B = 2.1266 + 0.66383s, 

D = 3.0301 - 1.7499s + 1.6466s2 , . , 

(A.20) 

E = 4.1282 + 1.6929s - 0.26292s^ 
E' = 0.89599 + 1.2761s - 0.15061s^ 
a = -0.78809 + 0.90278s. 

For the beauty-quark density 6had, for < 100 GeV^: 

a = -10.210 , /3 = -2.2296, 
A = -99.613 + 171.25s, B = 492.61 - 420.45s, 
D = 3.3917 + 0.084256s, E = 5.6829 - 0.23571s, 
E' = -2.0137 + 4.6955s, a = 0.82278 + 0.081818s. 

For the beauty-quark density 6had, for > 100 GeV^: 

a = 2.4198 , P = 0.40703, 

A = -2.1109 + 1.2711s, B = 9.0196 - 3.6082s, 

D = 3.6455 - 4.1353s + 2.3615s2, (A.22) 
E = 4.6196 + 2.4212s, E' = 0.66454 + 1.1109s, 
a = -0.98933 + 0.42366s + 0.15817s2. 
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FIG. 4: The same as in Fig. for > 20GeV2. 



31 





0.8 
0.6 
0.4 
0.2 


0.8 
0.6 
0.4 
0.2 



1 

0.8 
0.6 
0.4 
0.2 





L3 5.0 
TPC 5.1 ♦ 

TOPAZ 5.1 
DELPHI 5.2 -O 




AMY 6.8 
CELLO 7.0 ♦ 




01 



)PAL 8.9 
JPAL 9.0 ■♦■ 
PLUTO 9.2 y 



ALEPH 9.9 
OPAL 10.7 ■♦■ 
L3 10.8 




Delphi' 12.0 '-^ 

DELPHI 12.7 ■♦■ 



ALEPH 13.7 
OPAL 14.5 ♦ A P 



L3 15.3 
TOPAZ 16.0 ♦ A 




OPAL 17.5 
OPAL 17.8 ♦ 







0.5 



0.5 



0.5 



X 



FIG. 5: The same as in Fig. |21for a linear scale in x. 



32 



1 - 



0.5 h 



1.25 
1 

0.75 
0.5 
0.25 




_ 1 1 1 1 1 1 1 




~ CJKL 




; FFNS,,,, — - 




r GRS LO - - - - 




: SaSlD 





CELLh 26.8 




_L 



DELPhV 19.Q 
TPC 20.0 ♦ 
ALEPH 20.7 
DELPHI 21.0 <>- 




TASSO 23.0 
L3 23.1 ♦ 
JADE 24.0 



DELPHI 28.5 ^ 
OPAL 30.0 ♦ i 




DELPHI 40.0 
DELPHI 42.0 ♦ 



1.2 - 
0.9 '- 
0.6 
0.3 
^ 



PLUTO 45.0 ' - 

A i - 



ALEPH 56.5 - 
OPAL 59.0 ♦ 




AMY 73 
TOPAZ 80 ♦ 



1.5 
1 

0.5 


2 
1.5 
1 

0.5 






ALEPH 284 




0.5 



AMY 390 
DELPHI 400 ♦ 




DELPHI 700 
OPAL 780 ■♦ 



0.5 



0.5 



X 



FIG. 6: The same as in Fig. 0]for a linear scale in x. 



33 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 





1 — I 

CJKL - 
FFNSe,-M - 

GRS LO - - 
SaSlD - 
OPAL data ^ 




0.1 < rr < 0.6 



10 



100 



1000 



g2 [GeV^] 



FIG. 7: Comparison of the recent OPAL data for the Q^-dependence of the averaged^ over 
0.1 < x < 0.6 /a with the predictions of the CJKL and FFNScjfcZ models and GRS LO [^] and 
SaSlol^l parametrizations. 




1000 



Q2 [GeV^ 



FIG. 8: As in Fig. [7|for F2^(x,Q^)/a, averaged over four different x ranges. 



34 



0.5 

0.4 

0.3 

0.2 

0.1 


0.4 

0.3 - 

0.2 |i- 

0.1 


0.4 

0.3 

0.2 

0.1 





I ' ' ' I ' ' ' I ' 



g2 ^ 10 Gey2.-.-^-^^\ 




s 

GRV LO 
GRS LO 



SaSlD 

I . ■ . I . . . 



I ' ' ' I ' ' ' I ' ' ' I 

Q2 = 10 GeV2 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 


~ Q2 = 10 GeV^ ~ 




/ .* ~ 










* — ;r 




. . . 1 . . . 1 . . . 





0.2 0.4 0.6 0.8 1 
X 



0.5 
0.4 
0.3 
0.2 
0.1 




I ' ' ' I ' ' ' I ' ' 

xd^{x,Q^)loL 
g2 = 10 GeV2 




_L 



_L 



_L 



I ' ' ' I ' ' ' I 

xG^{x,Q^)loi 
g2 = 10 GeV2 




xh^{x,Q^)loL 
g2 = 10 GeV2 




0.1 

X 



0.2 



FIG. 9: Comparison of the parton densities predicted by various parametrizations at = 10 

GeV^, as a function of x. The charm- and bottom-quark distributions of the CJKL model are 

compared with the only other predictions for heavy-quark densities among the parametrizations 

used (given by the GRV LO and SaSlD parametrizations). 

35 



10^ 



102 - 




10^ r, 



10" - 



10- 



1000 GeV^ (x 100) 



CJKL 
FFNS,,-,; 

GRV LO 
GRS LO 



'--•..^^ 100 GeV2 (x 10) 




5 GeV2 (x 0.1) 




10-2 
10-^ 



10- 



10-2 

X 



10- 



10 



FIG. 10: Comparison of the up-quark distribution at four values of in the CJKL and FFNS 
models with the GRV LO J, GRS LO and SaSlD Q densities. 



36 




FIG. 11: The same as in Fig. EH for the gluon density. 
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FIG. 12: The same as in Fig. EH for the charm-quark density. 
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